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Abstract 

The Weierstrassian p, £ and a functions are generalized to R n . 
The n = 3 and n = 4 cases have already been used in instanton 
solutions which may be interpreted as explicit realizations of spacetime 
foam and the monopole condensate, respectively. The new functions 
satisfy higher dimensional versions of the periodicity properties and 
Legendre's relations obeyed by their familiar complex counterparts. 
For n = 4, the construction reproduces functions found earlier by 
Fueter using quaternionic methods. Integrating over lattice points 
along all directions but two, one recovers the original Weierstrassian 
elliptic functions. 
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Kaluza-Klein, Supergravity, String and Superstring theories and Super 
p-branes all involve dimensions beyond four, which then have to be compact- 
ified. The compactification results in lattices in the higher dimensions [|TJ. It 
might therefore be of interest to consider generalizations of doubly periodic 
functions to arbitrary R™. The R 4 case has already been treated by Fueter 
[|], who succeeded in obtaining quaternionic analogues of Weierstrassian el- 
liptic functions. In this note, we will present a unique and straightforward 
extension of Fueter 's results to R n . This general construction yields Fueter 's 
functions in n = 4 without requiring the use of quaternions. Furthermore, 
when one takes infinitesimal lattice spacings in all but two lattice directions 
and integrates over these points one recovers the usual Weierstrassian func- 
tions. 

Apart from possible future physics applications of n-tuply periodic func- 
tions, it is worth noting that (quasi)periodic classical solutions of Gravity and 
Yang-Mills theory have already been considered. Thus, inspired by Rossi's 
observation |§ that a singly- periodic instanton configuration represents a 
BPS monopole J|, Gursey and Tze || used Fueter's hyperelliptic functions 
to construct a solution with one Yang-Mills instanton per spacetime cell. This 
solution may be related to the monopole condensate || or to the Copenhagen 
vacuum 0, which is also based on a periodic arrangement of the theory's 
solitons. Hawking, on the other hand, argued M that metrics with a single 
gravitational instanton per unit cell make the dominant contribution to the 
path integral of Einstein's gravity. An adaptation of Weierstrassian a and 
( functions to three dimensions along the lines presented here leads to 
just such a self-dual solution, providing an explicit example of a "spacetime 
foam" based on Gibbons- Hawking multicenter metrics fTOf . 

We start by examining how the infinite sums defining the Weierstrassian 
functions are made to converge. These are 

CW = -/*,)* = I + £ { ^ + I + ±>, (2) 

2 

p{z) =\*<r{z) = I dz((z) = ln(z) + ^{ln( 2 -^)-ln(-^) + - + ^}, (3) 
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or, 

p(z) = \n(z) + £{ln(l "-) + - + f~ 2 } (4) 

In the above u = ri\uj\ + n 2 uj 2 , with (u^, uj 2 ) defining the basic lattice cell C. 
The parentheses {} ensure that each term in the sum is absolutely conver- 
gent; the series become meaningless if the parentheses are broken up and the 
terms they contain are separately summed. An efficient method for studying 
the convergence properties of (l)-(3) is provided by an integral test in which 
the double sum is replaced by S\£\ min d\uj\\uj\. This reveals why one, two and 
three subtraction terms are needed in (1), (2) and (3), respectively. Employ- 
ing a simple dimensional argument, one can proceed further and determine 
the precise forms of the subtraction terms, given only the first term in the 
sum (3), for example. To do this, we note (a) p is dimensionless (assigning 
the dimension of length to the coordinate z, say); (b) in the power series ex- 
pansion of (4), the highest power of uj permitted by convergence requirements 
is lu~ 3 . One then realizes that the subtraction terms following ln(l — ^) must 
be chosen so that p(z) has the expansion 

p(z) = \nz+Y,0(- 3 ) + ... (5) 

Hence using the power series 

ln(l - z/uj) = -z/uj - z 2 /2u 2 - z 3 /3uj 3 - .... (6) 

one chooses the three subtraction terms 

- ln(-cj) + z/uj + z 2 /2uj 2 (7) 

to arrive at the form (5). This will be the key in constructing R ra analogues 
of p(z). 

We will later need the well known facts that p(z) is doubly periodic while 
((z) and a(z) are quasiperiodic with the transformation properties 

C(z + u; h2 ) = C(*)+*7i,2 (8) 

and 

cr(z + ui t2 ) = -cr(z) expr]i )2 (z + lo 1)2 /2) (9) 
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where 771,2 = 2((cji i2 /2). 

We also record the so-called Legendre's relations 

I (( z )dz = 2ni (10) 

J dC 

and 

f ((z) dz = 2iri = r)icu 2 - r)2U 1 . (11) 

JdC 

We first rewrite (10) in the rather unconventional form reminiscent of 
Gauss's theorem 

JJ c dVV 2 p(z) = j> dB-Vp{z) = 27T, (12) 

where dV = dxdy and the "surface element" da = ndl. Here dl is the 
arclength and n a unit normal vector pointing outwards on dC. 

Let us next introduce the lattice basis vectors where a, p, = 1, ...,n. 
The volume of the unit cell C is given by 

V — —f fVi-VnJai) (a„) /io\ 

~~ n \ ai - an ■ y l6 > 

The basis vectors of the reciprocal lattice obey 

rfqf = r<°> ■ = 5 ab . (14) 
They are obtained from the q^ via 

r ("i) — _ r ^i-fn„(«2) „(«») (ie\ 
Vi (n-l)!\/ n " " ' 

We now seek higher dimensional versions of (12) in the form 

J dV n <9 M <9 M p n {x) = j> da^d^ p n (x) = - J du n , (16) 
where of course 

J du n = n n = 27r n/2 /r(n/2). (17) 

In order to obtain the result (16), the function p n (x), which is to serve as 
the analog of p(z), should have the form 

p n (x) oc G n {x) + X! •• X! G n (x — q) + (R n harmonics) . (18) 
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In (18), G n (x) is the Green's function for R n obeying 

d^d^G n (x) = -Q n S(x 1 )...S(x n ) = -VL n 5(x) (19) 

and 

q = n 1 q<- 1) + ... + n n q<- n \ (20) 

We have, of course 

GniX - q)= (2 O l — 2^ = Vx {H > 2 ^ (21) 

(ar — 2x • q + g 2 ) 2 |x — q\ 

The harmonics in (18) should now be chosen according to the general strat- 
egy outlined between equations (4) and (8). Thus in order to render con- 
vergent the sum (18), whose integral counterpart contains terms behaving 
like f\™\ min d\q\\q\ n ~ 1 /\q\ n ~' 2 i one again needs 3 subtraction terms. These are 
in fact nothing but the first three terms in the MacLauren expansion of (21) 
for \x\ « \q\min- This immediately leads to 

(22) 

Again, absolute convergence is attained only by considering each term de- 
fined by the outermost parentheses as an indivisible unit. Hence the terms 
proportional to q • x and to l/2g 2 cannot be summed separately (in which 
case they would appear to give zero by symmetry!) anymore than the uj^ 1 
term in (2) can. Note that the subtraction terms are indeed harmonics as 
anticipated in (18). 

This is perhaps an appropriate point for comparing our results with 
Fueter's. Fueter introduces the unit quaternions e M and corresponding 
to (I, ej) and (I, — e^), respectively. We have, as usual, e^ej = —5ij + e^k^u-i 
the indices i, j, k running from 1 to 3. Using these, one defines x = x^e^, x = 
x^e^jD = e M <9 M and D = e^d^. One then has x 2 = xx and DD = d^d^. 
Then, starting with the function 

n 1 V-r 1 1111111111, 

Z x = - + V{ + - + -x- + -x-x- + -x-x-x-}, 23 

x ^ x-q q qq qqq qqqq 

where q = n a q^ = n a q^e^, the q^ being lattice vectors, one arrives at the 
quaternionic version of Weierstrass's ((z) via 

C F (x)=^Z(x). (24) 
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We will see later that (24) indeed transforms the same way as ((z) under 
lattice shifts. The function corresponding to In a(z) is 

p(x) = £>Z(x) (25) 

and, remarkably, it turns out to have the quaternion-free form 

P {x) = - 2 + E(r^ -\-^r- ■ - ?V)}, (26) 

x q^o K x ~~ Q) Q Q Q 

coinciding with the n=4 case of (22). 

Returning next to (12), we see that the analog of ((z) is the n-gradient 

$\x) = dr Pn (x). (27) 

Taking n=4 and contracting both sides with the quaternion units e M , we 
recover the C F function of Fueter [2] 

C F = £V<9 M p 4 = £p4. (28) 

Leaving the quaternionic special case aside, let us now find the higher 
dimensional analogs of (8) and (9). Since 

d,d, p n = d, Ci n) = -n n £ .. £ S(x - q) (29) 

is a perfectly n-tuply periodic distribution, a shift by a lattice basis vector 
can only change by a constant vector. Thus 

C (n) {x + q (a) )=C (n) {x)+r] (n)(a } _ (3Q) 

This is indeed the transformation law for ( F when n=4. Putting x = —q^/2 
and noting that (j™\—x) = — (^(x), we find 

2Ci n W a) /2) = i n)(a) . (31) 

Integrating (30) and using p n { — x) = p n {x), we obtain 

Pn(x + q {a) ) = Pn (x) + ^ ■ (x + q^/2), (32) 

where there is no sum on the index (a) on the right hand side. Using (30), 
we can also obtain the higher dimensional form of Legendre's second relation 
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(11). Let us evaluate $ dC da^j^^x), where dC is the surface of the fundamen- 
tal hyperparallelepiped. Noting that the normal to the hyperplane defined 
by {q( a2 \ ..,q( an ^} is along the reciprocal vector given in (15), we find 

KE^ )(S) - W =^~ (33) 

a=l L \2) 

as the n-dimensional generalization of Legendre's second relation. 

There is an interesting relationship between p n and p n -\ which can be 
iterated until one finally arrives at the Weierstrassian elliptic functions. It 
is simplest to describe this relationship in the case of a rectangular lattice 
although it is generally valid. Starting with p n , one first chooses one lattice 
direction to coincide with, say, the x n axis and then lets the lattice spacing 
(in the same direction only) become infinitesimal. Then, integrating over q n , 
one obtains p n -i up to a multiplicative constant. Repeating this procedure 
until P3 is found, one finally has 

- r dq 3 p 3 = p(z)+p(z) = I dz((z) + / dz((z), (34) 

J — oo J J 

where the remaining lattice basis vectors (q^\q^) provide the periods 
{ui,U2). Thus the p n presented here are indeed intimately connected with 
Weierstrassian functions. We also observe that the split into an analytic and 
anti-analytic part is a feature of two dimensions, not shared for example by 
Fueter's quaternionic functions. 

Finally, the reader may wonder how the R n counterpart of p(z) is to be 
defined. It is clear that two derivatives of p n (x) will be involved. It is natural 
to classify the resulting functions according to their SO{n) transformation 
properties. One has the symmetric traceless traceless second rank SO(n) 
tensor 

= (d^d v - ^-dxdx) pn(x) (35) 

and the SO(n) scalar 

n (n) = dxdx Pn = -n n j2..T, 6 ( x - q) ( 36 ) 

already encountered in (29). Although both (35) and (36) are fully n-tuply 
periodic, the former involves true functions while the latter is an n-fold sum 
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over distributions, pffi is also similar to p(z) in that 

I da,p^= [ dV n d,p^ = (37) 



IdC 

just like p(z) which obeys 



p(z)dz = 0. (38) 
ac 
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